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We consider the fully overheated single-electron transistor, where the heat balance is determined 
entirely by electron transfers. We find three distinct transport regimes corresponding to cotunneling, 
single-electron tunneling, and a competition between the two. We find an anomalous sensitivity 
to temperature fluctuations at the crossover between the two latter regimes that manifests in an 
exceptionally large Fano factor of current noise. 



A single-electron transistor (SET) [1] , shown schemat- 
ically in Fig. 1(a), is one of the most thoroughly studied 
and widely used nanode vices. It has found its way to nu- 
merous applications in thermometry [2], single-electron 
pumping [3], charge detection [4, 5], and detection of na- 
noelectromechanical motion [6, 7]. The current noise in a 
SET limits the measurement sensitivity and is thus worth 
investigating [8]. 
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FIG. 1: (color online) (a) SET biased by a voltage Vbias- The 
charge in the central island can be tuned with the gate voltage 
Vg- (b) Coulomb diamonds in a symmetric SET. Blue dashed 
lines show the threshold voltage Vfa for SE tunneling, (c) 
Average temperature of the island versus bias voltage along 
the red vertical line in (b) for various tunnel conductances G. 
This illustrates the three transport regimes: I cotunneling, II 
competition, III SE tunneling. Dashed lines are asymptotes 
to pure cotunneling and pure SE tunneling. 



Nanodevices of sufficiently small size are overheated: 
the electron temperature in the device deviates from the 
lattice temperature. The temperature may fluctuate in 
this regime [9] , thereby affecting the current noise in the 
device. This motivates us to study overheating under 
Coulomb blockade conditions. In this Letter we concen- 
trate on a fully overheated SET. We assume the electron- 
phonon relaxation time, r e _ p ^, to exceed by far the elec- 
tron dwell time, so that the temperature is determined 
entirely from the balance of electronic heat transfers. In 
addition, we restrict our study to a symmetrically biased 
SET with junction conductance G <C Gq = e 2 /irh and a 
vanishing temperature of the leads. 

An early paper [10] addressed overheating in a SET 
in the regime of single-electron (SE) tunneling. It has 
been found that overheating instigates the SE transport 
at the threshold voltage Vc = (\/2— l)Vth, that is, well 
below the zero-temperature Coulomb blockade threshold 
V t h — 25E + /e, (dE* 1 — Eq T cVg is the charging energy, 
Eq = e 2 /2C) as shown in Fig. 1(b). We complement 
the consideration with electron cotunneling [11], which 
modifies the picture rather radically. We recognize that 
single-electron processes below V c try to cool the island, 
competing with the electron-hole excitations left behind 
by inelastic cotunneling that heat it up. This gives rise to 
a new transport regime: competition regime. The three 
regimes are evident in the voltage dependence of tem- 
perature as shown in Fig. 1(c). At low voltage cotun- 
neling dominates and the temperature scales with volt- 
age, fceT/eVb « 0.24, as expected for a fully overheated 
nanode vice. Sufficiently high temperature activates SE 
transfers that cool the island and sets the competition 
in. The temperature/voltage ratio reaches the minimum 
k B T c /eV c = l/V21n(const. x G Q /G) near V c . Above 
the threshold, the SE processes heat the island result- 
ing in an increase of temperature. A pure SE picture 
captures only this rise predicting T = Q at Vt> < Vc- 

The most interesting features can be found at Vy, « Vc 
where the crossover between competition and SE regimes 
takes place. We show that near the crossover the electric 
current is anomalously sensitive to temperature changes: 
it is significantly modified by a temperature change 5T <C 
T. The underlying mechanism of the anomalous sensitiv- 



ity is the strong temperature dependence of thermally ac- 
tivated tunneling rates. The overheated SET also detects 
fluctuations of its own temperature, manifest in an en- 
hanced current noise. The current noise Sj is commonly 
characterized by Fano factor F = Si/2el, F < 1 for 
most of nanode vices. At the crossover, F oc ln 4 (G Q /G) 
reaching an impressive 1500 for G/Gq — 10~ 6 . 

We implement a method that allows to access the full 
statistics of temperature and current fluctuations. The 
statistics are described with an action S depending on 
counting fields x, £, conjugated to the transferred charge 
and the energy of the island, respectively [9, 12]. For 
single-electron tunneling, the dynamics of a SET is gov- 
erned by a master equation: The stationary probability 
distribution of charge states labeled by N, p^, satisfies 
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where elements of Sa^at' correspond to single-electron 
tunneling rates so that N' — N ± 1. It is shown in the 
theory of full counting statistics (FCS) [12, 13] that in 
order to obtain the action, one should modify J^n,N' to 
include counting fields \y £■ The action is then given by 
the eigenvalue of so-modified matrix with the smallest 
real part, <t(x,£) [13]. One could include higher-order 
tunneling processes by replacing Hn,n' with self-energies 
composed of all possible irreducible tunneling diagrams 
that take the SET from charge state N' to N [14]. This 
is the way to account for cotunneling. Usually, if the 
tunneling processes of different orders become equally 
important, the situation is very difficult to comprehend 
[15]. This situation typically occurs if the rate of electron 
transfer is comparable with the energy released in the 
course of transfer. This implies that the flow of charges 
can not be divided into separate events of any order. 

Fortunately, this is not the case of the fully overheated 
SET where cotunneling and SE events are separated even 
if they are equally important for transport. To under- 
stand this, let us concentrate on a blockaded diamond 
corresponding to a certain charge state, say N = 0, and 
assume T <C eVth/ks- The first SE transfer must be ther- 
mally activated and proceeds with the suppressed rate 
rexp(-W7fc B T), W = e(V th -V b )/2. It brings the island 
to the closest excited state N = 1. However, the island 
will quickly get back to N = 0: The first SE transfer 
is followed by a second, after a time given by an unsup- 
pressed rate T. Similarly, a cotunneling event can also 
be viewed as two SE events separated by a time interval 
fi/ eVth [15]. We see that the transport separates to ele- 
mentary events each encompassing two SE transfers. The 
events are independent since time interval between them 
exceeds the time separation between the transfers by a 
large factor, min(Gq/G, exp(W/k B T)). Therefore, the 
cotunneling and SE contributions can simply be summed 
together, S = S sc + S cot . 
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FIG. 2: (color online) (a) Differential conductance as a func- 
tion of Vb for several values of G. Dashed line is an asymp- 
tote for pure SE tunneling. (b) Temperature sensitivity 
d(ln/)/d(lnT) as a function of Vb for the same values. Inset: 
Fraction of SE transfers in the current flow. 



Analytical results can be obtained by taking into ac- 
count only two charge states on the island, N — and 
N = 1. However, the validity of this approach requires 
log(Gq/G) ^> 1, rarely the case in practical devices. 
Therefore we also perform accurate numerics, where we 
take more charge states for S sc and weight 6> co t with the 
probabilities of those states. 

Let us first outline the three different regimes men- 
tioned. To simplify the formulas, we set e = k-Q = h = 1 
and define 7 = G/ttGq. The SE part of the action in the 
relevant limit T <C Vb, V"o; T£ <C 1 reads (see Appendix) 
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Using this action one evaluates the charge current, 
I = d x S\ x =n, and heat current, H = d^S\^ = Q, as 
functions of temperature. Equating the latter to zero 
yields the average temperature in the SE regime, T — 
4^ (Vb - V c ) (V b + V c + 2V th ) above V c . The current 
steeply rises at the threshold, I sc oc exp(— Wc/T) — 
exp(-^) at V b «V c . 

The cotunneling regime takes place at Vb <C Vq where 
T w Vb. In this region, the action reads (see Appendix) 



S co t =7 



ttT 3 [ ( 



1 1 

SE+ + SE 
+2J(0, x) - X(v, 0) - X{-v, 0) - 2J(0, 0)] 



x X(v, x) + e x I(—v, x) 
(3) 



where l(v,x) = dz 4(z _ io e +y sinh . (<+Ue) , v = Vb/vrT, 
and x = 7rT£. This yields the average temperature T w 
0.238Vb and electric current 
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Let us note that the current in the absence of overheating 
is also oc and is given by the same expression with 
coefficient 1/12-7T. Therefore, the overheating enhances 
cotunneling current roughly by a factor of 2. 

Thereby we resolve a long-standing discrepancy be- 
tween theory and experiment. The pioneering work [16] 
on cotunneling in SET has reported such a factor of 
2 mismatch for the most conductive junctions. This 
is explained by full overheating. For the less conduc- 
tive junctions the mismatch factors were 1.4 and 1.2, 
explained by incomplete overheating. In this case the 
smaller electronic heat flows may have become compa- 
rable with phonon heat transfer. In addition, Rcf. 16 
reports a crossover to SE regime at approximately half 
of the Vth expected: this conforms to theoretical value of 
V c . 

A further increase of Vb increases the temperature 
and activates SE processes. Comparing cotunneling 
(~ "fV£/W 2 ) and SE (~ jT exp(-W/T)) rates, we ex- 
pect the SE processes to become important at V\, ~ T ~ 
IVln(l/ 7 ) <C Vc- We enter the regime where equilib- 
rium temperature is determined from the competition 
between SE transfers cooling the island and cotunneling 
events heating the island. Similar estimation gives with 
logarithmic accuracy T w W/hi(l/7) in the whole in- 
terval of the competition, that is, up to Vc- Inset in 
Fig. 2(b) shows the fraction of SE events in the current 
flow. The fraction grows almost linearly from the border 
of the cotunneling regime up to Vc- Indeed, at low volt- 
ages each cotunneling process provides Vb of heat while 
a SE transfer cools the island by a value of ~ W: Many 
cotunneling events match up a single SE transfer. Near 
Vc, a SE transfer gives a vanishing cooling ~ (Vc — Vb): 
Many SE transfers are needed to balance a single co- 
tunneling. A simple analytical expression for the total 
current in this regime does not exist. Qualitatively, it is 
estimated by I = (1 + (1- Vb/Vc) _1 )/cot, ^cot oc 7 2 being 
the cotunneling current in the absence of overheating. 

It is interesting to note an anomalously high tem- 
perature sensitivity of the fully overheated SET. The 
underlying mechanism is the same as in living organ- 
isms that rely on the balance of the thermally acti- 
vated rates ~ exp(— W/T). Changing a rate by a fac- 
tor of e is achieved by a small temperature change 
ST ~ T(T/W) < T and may even lead to the destruc- 
tion of an organism. Similarly, the thermally-activated 
character of SE transfers may result in a high sensitiv- 
ity of the current that we characterize by a dimension- 
less number d(ln/)/d(lnT) ~ W/T. This is plotted in 
Fig. 2(b). We see that the sensitivity reaches the maxi- 
mum at the crossover between competition regime where 
it indeed scales as Wc/Tc- The linear growth below Vc 
is explained by the almost linear increase of the frac- 
tion of SE electron transfers in the competition regime. 
Above Vc, the sensitivity drops like l/(Vb — Vc) owing 
to temperature increase. 



Discreteness of charge transfer through the structure 
gives rise to a current noise, Si, and heat current noise, 
Sfj, both white at frequencies u -C /. The heat cur- 
rent noise produces temperature fluctuations that persist 
over a significant time, r = (dJ r /dT)/(dH/dT) (T be- 
ing the total free energy of the island, proportional to 
its volume). The fluctuation of temperature is given by 
((AT) 2 ) = S h /T(dH/dT) 2 . Temperature fluctuations 
change the current, 51 = (dI/dT)5T(t) giving rise to ex- 
tra "slow" current noise persisting at frequencies ~ 1/t, 

Si - s1ow = (Sylr) s »> (4) 

manifesting the overheating. Anomalous temperature 
sensitivity gives rise to an anomalous Fano factor, plotted 
in Fig. 3. Similar to sensitivity, the Fano factor also peaks 
at the crossover between SE and competition regimes. 

Let us now concentrate on the crossover region. Three 
factors contribute to the heat balance at Vb ~ Vc: 

= H cot + 2- 1 l 2 1 Te- w ' T (V h - V C ) 7 T 2 e~ w / T , (5) 

cotunneling heating H cot oc 7 2 V^ that is approximately 
constant, SE flow that switches from cooling to heating 
at Vb = Vc, and extra cooling that stabilizes tempera- 
ture in SE regime. We define the temperature Tc at the 
crossover through 

gcotVc / Vc _M =1 .Ic 2- 1/2 
V2~fT c ° XP \V2Tc V2j V C ~ In(l/7) ' 

and introduce dimensionless deviations of voltage x 
and temperature y such that Vb = Vc + V%~Tc + 
2xVc(T c /Vc) 2 , T - T c + \/2yVc(Tc/Vc) 2 , (T c /V c ) be- 
ing an important dimensionless small parameter enabling 
the scaling. The heat balance rescales to 

e y (x - y) + 1 = 0, (6) 

which implicitly gives the temperature as a function of 
voltage. The crossover takes place at x, y ~ 1, the 
rescaled equation being valid in a larger interval up to 
x,y s=s Vc/Tc- We see that the crossover is shifted from 
V c by V2T C - The width of crossover interval in volt- 
age/temperature is small ~ T c (Vc/T c ) <C T c . While 
temperature changes insignificantly, the relative change 
of quantities of interest is by an order of magnitude. Cur- 
rent is rescaled to 

Similarly, 

— = V2(1 + e-y)I(y). (8) 
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FIG. 3: (color online) Fano factor of the temperature fluc- 
tuation induced current noise for different values of G/Gq. 
Dashed line is an analytic approximation for G/Gq = 10~ 6 , 
and agrees well with the numerical result. Inset shows the 
Fano factor for G/Gq = 10 -6 in an absolute voltage scale for 
various values of eVa/Ec- The peaks fall on top of each other 
once rescaled to common Vc 



This yields the Fano factor at the crossover 



p {dl/dTf S ti 
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Vc 



(1 



(9) 



describing a sharp rise as y — > 0. Its fall in the SE region 
is described by substituting Tc — > T(Vb). This fits well 
to the numerical results as shown in Fig. 3. 

There are very interesting statistics of temperature 
fluctuations at the crossover. Generally [9], one expects 
deviations from Gaussian statistics for temperature de- 
viations of the order of average temperature that oc- 
cur with exponentially small probability InP ~ —T/8s, 
5s being the single-electron level spacing in the island. 
In an overheated SET around the crossover, the de- 
viations 5T ~ Tq(Tc/Vc) are already non-Gaussian 
and their probability is greatly enhanced, In P ~ 
-(Tc/Ss)(T c /Vc) i - We will present the detailed results 
for the distribution of the fluctuations in a separate publi- 
cation. The system ideally suits for the experimental ob- 
servation of temperature fluctuations at nanoscale: large 
fluctuations are easily read by a measurement of the elec- 
tric current, slowly fluctuating in time. 

Finally, let us estimate the importance of electron- 
phonon interaction to assess the feasibility of full over- 
heating. For such an estimate, it is enough to add a 
term to the action S e - ph = -£VT 5 £(1 - 5Tf), where V 
is the volume of the island and S the material-specific 
electron-phonon coupling constant. The island is over- 



heated provided SVT 3 <C 7- For typical values [17], £ ss 
2 x lCT 9 WK-'Vm -3 , T w 0.1 K, and 7 w 1(T 3 , the vol- 
ume of the island should be of the order of V « 10 -4 /jm 3 
to reach this regime. This is easily achievable experi- 
mentally. High Fano factor sw 10 3 requires 7 w 10 -5 
and V ~ 10 -6 /im 3 , feasible in smaller systems such as 
granular metals and multi- walled carbon nanotubes. 

To conclude, we have studied the fully overheated SET 
revealing the importance of cotunneling processes that 
compete with single-electron transfers in a wide inter- 
val of bias voltages. The fully overheated SET exhibits 
anomalous temperature sensitivity and "slow" current 
noise with a huge Fano factor as a result of temperature 
fluctuations. These effects are most pronounced at the 
crossover between competition and single-electron tun- 
neling dominated regimes. 
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DERIVATION OF EQ. (2) 



Let us define a function 



T -A(|+l/T) 

T a ( x , A) = 7"e-* 1 + 2^(1, 1 + 2 + TO - e - A/T ), 



(10) 



where 2-Fi(a, b; c; z) is the Gauss hypergeometric function. The counting field modified single-electron tunneling rates 
at the left junction are 

f N +1 , N = f l ( X l - XD,8E+ -L1 L + n D ), t L N _^ N = f L ( X D - XL,SE N + hl- Ho), 

and similarly for the right junction. When T <C |A| we have 



2 Fi(i, 1 + to 2 + to - e - A/T ) « j gA/T 



1 _ i±n a/t 

1 2+T{ e ' 

e^r(-ror(2 + to + "^g^ 



A > 
A < 



f Q (x,A) 



7«e-^p^, A>0 

7 a e - X e A *-l ; A< q 



(11) 
(12) 



Here T(z) is the gamma function. 

At low temperatures it is enough to consider two charge states, N = and N = 1, on the island. The effective 
action is then easy to calculate, yielding 



r ,i + r li0 - yj (r ,i + i\ ) 2 + 4(f ,ifi, - r ,iri, ) 



rcuiYo - ro,iri ; o 

To.i + Tin 



(13) 



The last approximation is valid when r ,i ^> I\ or r ,i <C I\ , i.e., under Coulomb blockade. Using the form of 
Eq. (12) this becomes (for xl = X, Xd = Xfl = 0, 7 L = 1 R = 7, Ml = = eV h /2, and 5E+ = E c - eV G ) 



S = - 



jTe -SE+/T 



e(i + n)SE+ 1 + j|+ e -^+/T C0Sh 



cosh 



2 V T 



2£ -X 



-e- 5B+ «cosh 



eUb /l 



e-^ + «cosh 



eU b 



2 V T 



2 V T 

This results in Eq. (2) of the main text in the limit T <C V b , V c ; T£ < 1. 



I ' 27 ' ) [1-2^(1 + TO] 



(14) 



DERIVATION OF EQ. (3) 

Cotunneling rate with the counting fields between leads a and [3 (<G {L, R}) held at zero temperature 



1 



+ eH/2 - <5T+ -E - eF b /2 - 5E N , 

/0 /'00 
de[l - f(e - E)}e^ e - E) / de'/(e' — E — eVb)e~^ e ' ~ E ~ eVb \ (15) 
-oo Jo 

where /(-E 1 ) is the Fermi function at the island temperature and eVb = /U a — [ip. Now write 

/oo />oo 
dtfi^e 1 ^-^*, f{e' -E-eV h )e-^'- E - eV ^ = dt 1 F 2 {t')e^' - E ~ eV ^' , (16) 
-oo J —00 



and assume eVf, <C <5-E^, i^iu so that 



cot ' N ~ 7 7 U# + «5j loo (t-iO) 2 



(17) 



G 



Carrying out the Fourier transforms gives us 



Fl{t) 2sinh(7rT(i + i£))' ^ ~ 2 sinh(7rT(t - i£)) ' 



resulting in 



2 „oo it)(z-iO+) 

r »-=^ V fe + ^J jT3 L d - 4(,-iO + )^^( g + i x) - (19) 

where v = eV^/nT and x — irTt;. This leads to Eq. (3) of the main text when the four different processes, corresponding 
to cotunncling from the left lead to the right lead, right to left, left to left, and right to right, are summed and the 
X = £ = part is subtracted. The integral above can be evaluated by closing the contour around the upper (lower) 
half plane for v > (v < 0). The result is 

. 7rw + 2cota; tt ^ e vx ~ nnv (vx - nirv - 2) 
2 sin 2 a; 2 ^ (x - mr) 3 

„ . 7T ^ e ra+m ™(TO + 717TV - 2) 

J(«<0,ar) = -- V ^ -. 21 

v ; 2 ^ (a; + n?r) 3 v ' 

n=l K ' 



These sums can be used to evaluate the action to a desired accuracy. 



